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ABSTRACT 

We show that changes in the orientation of the inner accretion disk of Cygnus X-l affect the shape of the broad 
Fe Ka emission line emitted from this object, in such a way that eV-level spectral resolution observations (such as 
those that will be carried out by the ASTRO-E2 satellite) can be used to analyze the dynamics of the disk. We here 
present a potential diagnostic tool, supported by numerical simulations, by which a few observations of Cygnus 
X-l, separated in time, can determine whether its accretion disk actually precesses, and if so, determine its period 
and precession angle. This approach could also be used for similar studies in other microquasar systems. 

Subject headings: X-rays: binaries—X-rays: individual (Cygnus X-l) 


1. INTRODUCTION 

Cygnus X-l is a 5.6-day X-ray binary that harbors the best 
studied black-hole candidate in the Galaxy. The mass of the 
accreting compact object has been estimated as ~ 10.1 M 0 and 
the donor star is classified as an 09.7 lab supergiant of ~ 17.8 
Mq (Herrero et al. 1995). The system is located at ~ 2 kpc 
(e.g. Gierlinski et al. 1999). The black hole accretes through 
the wind of the companion star. Most of the time, the X-ray 
source is in the so-called low/hard state, characterized by a rel¬ 
atively weak blackbody component peaking at a few keV plus 
a strong hard power-law of photon index ~ 1.6. A nonthermal 
radio jet has been observed in this state (Stirling et al. 2001), 
extending up to ~ 15 mas. The jet seems to form an average 
angle with the line of sight of ~ 30° (Fender 2001) and it has 
been suggested that it might be precessing (Stirling et al. 2001, 
Romero et al. 2002). Occasionally, a transition to a high/soft 
state can occur. In this state, most of the radiated energy is 
concentrated in the blackbody, while the power-law component 
becomes softer, with an index of ~ 2.8 and no jet has been ob¬ 
served. 

The usual interpretation of the X-ray behavior of the source 
is that the blackbody component originates in a cold, optically 
thick accretion disk, whereas the power-law component is pro¬ 
duced in an optically thin hot corona by thermal Comptoniza- 
tion of disk photons (Poutanen et al. 1997, Dove et al. 1997, 
Esin et al. 1997, 1998). The hot corona fills the inner few 
tens of gravitational radii around the compact object and the 
accretion disk would penetrate only marginally into the coronal 
region. In the low/hard state the thermal X-ray luminosity is 
dominated by the corona, with typical luminosities of ~ a few 
times 10 37 erg s -1 . During the transition to the high/soft state 
the corona is likely ejected as the accretion disk approaches to 
the black hole (Fender et al. 2004, also Esin et al. 1997, 1998). 
Most of the energy is then dissipated by the disk, until the in¬ 
ner part of it dominates the radiation again, and the cycle starts 
again. 

In the low/hard state, the disk is illuminated by hard photons 
from the corona resulting in the production of an Fe Ka line 


and a Compton reflection feature. The first detection of the line 
was made by Barr et al. (1985) with EXOSAT. They reported 
a broad (FWHM ~ 1.2 keV) emission line at ~ 6.2 keV with 
an equivalent width of ~ 120 eV. Kitamoto et al. (1990) ob¬ 
tained a Tenma GSPC spectrum which was consistent with a 
narrow emission line at ~ 6.5 keV with an equivalent width of 
60-80 eV. A Compton reflection feature was then found above 
20 keV (see Tanaka 1991 and references therein). The Ginga 
spectrum in the 2-30 keV range can be fitted quite well by the 
sum of a power-law with index ~ 1.7, a reflection component 
and a narrow Fe emission line at 6.4 keV with an equivalent 
width of ~ 60 eV (Tanaka 1991). Subsequent ASCA observa¬ 
tions confirmed these results but restricting the width to 10-30 
eV (Ebisawa et al. 1996). A broad edge at E > 7 keV was also 
reported. A detailed historical account of attempts to detect an 
Fe line in Cygnus X-l and the intrinsic difficulties it entails, can 
be found in Reynolds & Nowak (2003). Recently, Chandra ob¬ 
served Cygnus X-l with the High Energy Transmission Grating 
Spectrometer in an intermediate X-ray state (Miller et al. 2002). 
The narrow Fe line was detected at E = 6.415 ±0.007 keV with 
an equivalent width of W = 16^ 3 eV, along with a broad line 
at E = 5.82 ±0.07 keV with W = 140^ eV. A smeared edge 
was also detected at 7.3 ±0.2 keV. Miller et al. (2002) interpret 
these results in terms of an accretion disk with irradiation of the 
inner disk producing the broad Fe Ka emission line and edge, 
and irradiation of the outer disk producing the narrow line. The 
broad line is shaped by Doppler and gravitational effects and, 
to a lesser extent, by Compton reflection. 

For different spectral states, the different disk structure may 
change the Fe Ka line. For example, if the disk is truncated 
at a much larger radius rather than the innermost stable orbit 
as suggested for the hard state, then the width of the Fe Ka 
line may significantly decrease. In any case, the variation due 
to disk precession, if that were observable, and that due to the 
accretion mode (disk structure), will have a different temporal 
signature (periodicity), what would make them easy to distin¬ 
guish. 

In this paper, we show that changes in the orientation of the 
inner accretion disk of Cygnus X-l would affect the shape of 
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the broad Fe Ko emission line in a periodic way. Under the 
assumption that the X-ray spectrum of Cyg X-l is not substan¬ 
tially more complex than what has already been found in the 
Chandra observation, eV-level spectral resolution observations 
of the system (as those that will be carried out by the ASTRO- 
E2 satellite) can be used to constrain the dynamics of the disk 
and establish whether the accretion disk of Cygnus X-l actually 
precesses, and if so, to determine its dynamics. 


2. DISK PRECESSION IN CYGNUS X-l 


Several X-ray binaries present periodic behavior in their light 
curves on timescales longer than the orbital period. Among 
these systems we can mention Her X-l, SS 433, and LMC X- 
4. It has been suggested that these long periods correspond to 
the precession of the accretion disk (e.g. Katz 1973). In the 
case of SS 433 the precession is directly measured in the jets, 
so if these are attached to the accretion disk it is reasonable to 
expect that the disk will also display precession (Katz 1980). 
Although there is no reported compelling evidence yet for disk 
precession in black hole binaries, it is reasonable that the same 
mechanisms responsible for this phenomenon in neutron bina¬ 
ries will apply. 

The mechanism that produces the precession might be the 
instability of the response of the disc to the radiation reaction 
force from the illumination by the central source (e.g, Wijers 
& Pringle 1999, Ogilvie & Dubus 2001), or the tidal force of 
the companion star on a disk which is not coplanar with the 
binary orbit (Katz 1973, Larwood 1998, Kaufman Bernado et 
al. 2002). Uniform disk precession will occur in this case 
only if the sound crossing time through the disk is consider¬ 
ably shorter than the characteristic precession period induced 
by the perturbing star. The precession angular velocity is given 
by (e.g. Romero et al. 2000): |fi p | ss |p^^cos0, where G is 
the gravitational constant, r m is the orbital radius, Wd is the inner 
disk angular velocity, 9 is the half-opening angle of the preces¬ 
sion cone, and m is the mass of the star that exerts the torque 
upon the disk. The orbital period T m is related with the involved 

masses and the size of the orbit by Kepler’s law: r^ n = , 

where M is the mass of the accreting object. The ratio between 
the orbital and the precessing periods can be related through the 
disk angular velocity wa = (GM/fj) 1 / 2 : 
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where T v = 2it/Cl v and k = ra/r m . Since k < 1, normally 
T m /Tp < 1. In the case of Cygnus X-l, Brocksopp et al. (1999) 
have reported multiwavelength evidence for the existence of a 
period of 142.0 ±7.1 days. Similar precessing periods have 
been calculated by Larwood (1998), Katz (1973, 1980) and 
Katz et al. (1982) for other X-ray binaries on the basis of the 
same model. As shown by Romero et al. (2002) in the case 
of Cygnus X-1, such a period can be obtained from tidally- 
induced precession for an accretion disk with a size ~4x 10 11 
cm, if the half opening angle of the precession cone is ~ 15°. 
For a purely wind-fed system, this size might be too large and 
other mechanism may be in operation to generate the observed 
timescales. In particular, radiation-driven precession (Pringle 
1996, Maloney & Begelman 1997, Ogilvie & Dubus 2001), 
wind-driven warping and precession (Schandl & Meyer 1994, 
Quillen 2001), and spin-spin precession (Bardeen & Petterson 
1975, Armitage & Natarajan 1999) can yield precession peri¬ 
ods of several weeks to a few months. In all these mechanisms 


the observed optical and X-ray modulation points to precession 
of the disk, whereas the radio variations might originate in the 
jet. 


3. FE K-a LINE PROFILE DIAGNOSIS 

The use of emission lines as a diagnosis tool for the state of 
binary systems has been proposed in the past (e.g., for an in¬ 
vestigation on supermassive black hole binarity, see Torres et 
al. 2003, also Gaskell 2003; Zakharov et al. 2004a,b). Here we 
show that a similar method can be used to extract information 
about the precession status of microquasars. 

For the case of Cygnus X-l, we shall assume that the time- 
averaged disk inclination angle is 35°, which is in agreement 
with the fitting of the system’s Fe line (Miller et al. 2002). 
Values around this time-averaged inclination angle were also 
found for other binary systems ( e.g.. Her X-1, LMC X-4, SMC 
X-l, etc. e.g., as discussed by Larwood 1998). We also as¬ 
sume two extreme cases for the amplitude of the disk preces¬ 
sion: (1) the disk inclination angle precesses from 31° to 39° 
(very low magnitude of the precession angle); (2) the disk in¬ 
clination angle precesses from 5° to 65° (large magnitude of 
the precession angle). The amplitude of the precession of the 
inner disk should not be large if it is due to the tidal force of 
the secondary star and if the disk (especially the outer disk re¬ 
gion) develops a significant warp. However, if the initial spin 
direction of the BH is significantly different from the orbital 
angular momentum direction, the inner disk, which is confined 
to the equatorial plane of the BH (if the spin is high) due to the 
Bardeen-Peterson effect, may precess around the total angular 
momentum (dominated by the orbital angular momentum) with 
an amplitude as large as the initial orbital inclination angle with 
respect to the BH equatorial plane. If we assume that the disk 
is rigidly precessing around the total angular momentum (dom¬ 
inated by the orbital angular momentum), then the amplitude of 
the precession may also be around 30°. 

Several calculations on the disk line profiles have been per¬ 
formed. We use a ray-tracing technique and elliptic integrals 
(Rauch & Blandford 1994; see also Yu & Lu 2000; Lu & Yu 
2001 and references therein) to follow the trajectories of pho¬ 
tons from the observer, keeping track of all coordinates until 
the photons either intersect the accretion disk plane, disappear 
below the event horizon, or escape to "infinity" (operationally 
defined to be r = 1000 GM/c 2 away from the BH). We then cal¬ 
culate the redshift factor for a photon (to the observer) emitted 
from a particular position on the disk. The solid angle sub¬ 
tended at the observer by each disk element is also calculated. 
We set the inner radius of the disk to be at the marginally sta¬ 
ble orbit (6r g for a Schwarzschild black hole or 1.23r g for a 
Kerr black hole with spin a/M = 0.998, where r g = GM/c 2 ), 
and the outer radius at 160r g . We assume that the surface emis- 
sivity of line photons follows a power-law, r~ q , with q = 2.5. 
Both the power-law emissivity law and the size of the disk in 
Schwarszchild units, are usual assumptions (see, e.g., Nandra 
et al. 1997). The BH spin is assumed to be a/M = 0.998. In mi¬ 
croquasar systems, both the high frequency quasi-periodic os¬ 
cillation and relativistic lines suggest a high spin. In any case, 
we proved that if we were to assume a lower spin, there is not 
much qualitative difference for the problem we have studied 
here. With the above assumptions, we sum up all the photons 
received by the observer, which is emitted from each disk ele¬ 
ment, and obtain the profile of emergent Fe Ka lines, with dif¬ 
ferent inclination angles from the Cygnus X-l system, as shown 
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in the Figure 1. 

4. OBSERVING THE FE LINE PROFILE VARIATIONS 

The complex X-ray spectrum of Cygnus X-l revealed by 
the Chandra observation reported in Miller et al (2002) im¬ 
plies that the detection of minute variations in the Fe line shape 
will require a high spectral resolution instrument with large 
throughput. In what follows, we discuss the feasibility of de¬ 
tecting the precession of the accretion disk with ASTRO-E2, 
and specifically with the X-ray Spectrometer (XRS) consisting 
of an array of semiconductor-based calorimeters delivering the 
best spectral resolution to date at 6 keV (pre-flight value of 6.5 
eV). ASTRO-E2 is a JAXA/NASA mission to observe X-rays 
with unprecedented high spectral resolution imaging detectors, 
which is scheduled for launch by Summer 2005. The ASTRO- 
E2 Science Working Group target list includes two observations 
of Cygnus X-l. 

The underlying assumption is that the X-ray spectrum of Cyg 
X-l can be well described by the best-fit model obtained by 
Miller et al (2002). Indeed, at higher spectral resolution, other - 
so far undetected- spectral components might appear. Examples 
include the presence of blue-shifted edges due to high velocity 
ejecta. Should these, or other, features impact strongly either on 
the Fe line region itself or in the energy ranges where the un¬ 
derlying continuum is estimated, the specific methodology pro¬ 
posed here will need to be revised. Fortunately the instrument 
selected will deliver spectra with a very high signal to noise at a 
spectral resolution so high that many of these putative compo¬ 
nents might be resolved and modeled out. That would certainly 
complicate, but not invalidate, the proposed analysis. 

As argued by Miller et al (2002) the spectrum of Cygnus X-1 
is complex below 3 keV, where they did not succeed in fitting 
an appropriate model. All our discussion assumes the Be filter 
is on at the XRS. This removes low energy photons which are 
not needed for our purposes. In our analysis, we ignore photons 
below 3 keV and above 9 keV. In addition, we do not use the X- 
ray Imaging Spectrometers (XIS) on board ASTRO-E2 , which 
have less spectral resolution at the Fe line energy range. A typ¬ 
ical 50 ks exposure is assumed, and the background is assumed 
to be negligible for such an strong source (see below) and not 
included in the simulations. Pre-launch calibration redistribu¬ 
tion matrices and efficiency curves have been downloaded from 
the ASTRO-E2 web site at NASA 1 . 

The simulated model is that fitted by Miller et al. (2002) to 
the Chandra data. The continuum is a power law with T = 1.8, 
absorbed by a Galactic column of Nh = 8.1 x 10 21 cm -2 which 
is absorbed by a smeared edge at 7.2 keV, with a width of 7 
keV and a depth of 1.2. The narrow line component at 6.415 
keV, believed to arise from the irradiated outer disk, has been 
simulated as a gaussian of width tTnarrow = 30 eV. The broad line 
component has been simulated using the numerical models ex¬ 
plained in the previous section, for a fixed equivalent width of 
140 eV and a variety of disk inclination angles. 

As expected, this model produces a very high count rate in 
the XRS, of the order of 60 counts s _1 . This count rate will 
be distributed among several XRS pixels, according to the PSF. 
Although the overall count rate is below the telemetry limit, the 
fraction of events that will be measured by the on-board soft¬ 


ware as medium or low resolution will be large (~ 30-40%). 
There are two possibilities to deal with this: either using the 
neutral density filter (which will decrease the overall count rate 
to an acceptable level of ~ 6 counts s _1 ) or to ignore the few 
pixels with the higher count rates and work only with the pixels 
which have count rates below a few counts s _1 . Both of these 
procedures will result in an overall loss of throughput. This is 
why we consider effective exposures of 10 and 5 ks. 

To analyze the simulated data, we follow Miller et al (2002) 
to fit the continuum, by excluding the range from 4.0 to 7.2 ke V. 
A single power law leaves enormous residuals which can be 
well fitted by the smeared absorption edge. In general, the edge 
energy is well reproduced by the fit (statistical 90% errors in the 
range of 50-150 eV depending on the effective exposure time), 
although there is substantial degeneracy between the width of 
the edge and its depth. This does not affect the continuum in 
the Fe emission line region. 

Once the continuum is fitted, we include all data in the 3.0- 
9.0 energy band, and add a narrow gaussian emission line and 
a relativistic emission line (see Figure 2, left). Thanks to the 
superb spectral resolution of this instrument, the narrow line is 
very well characterized, with errors in its centroid actually lim¬ 
ited by systematics (2 eV) rather than by statistics, even in a 5 
ks exposure. The relativistic line model returns the disk incli¬ 
nation angle with a 90% error of 0.3 deg for a 50 ks exposure 
and 0.7 deg for a 5 ks exposure. This is due to the fact that the 
sharp drop in the blue edge of the line is very clearly marked by 
the XRS. Figure 2 (right panel) shows the differences in these 
sharp edges for a 50 ks exposure and various disk inclination 
angles. Note that these changes amount to about 50 eV per de¬ 
gree of disk inclination, and therefore the claimed limit in the 
systematics for line centering of 2 eV is really not an issue for 
this purpose. Table 1 summarizes the results of the fits to simu¬ 
lated data with 5, 10 and 50 ks net exposure and disk inclination 
of 35°. 

5. CONCLUDING REMARKS 

This work suggests a diagnostic tool to investigate whether 
the accretion disk of Cygnus X-l (and of other microquasar 
systems) precesses, and if so, what is the period and precession 
angle. We have shown that the study of the periodic variations 
of the Fe Ko line, that would be unavoidably produced in the 
putatively precessing disk of the system, are observable in short 
(5-10 ks) exposures of the ASTRO-E2 satellite. The degree of 
precision and confidence level up to which we will be able to 
determine the inclination angle of disk for each short obser¬ 
vations, thus, the magnitude of the precession, was shown to 
be sufficiently high as to allow a clear determination of these 
parameters even when the precession angle is of only a few de¬ 
grees. 

The work of DFT was performed under the auspices of the 
U.S. DOE-NNSAby U. of California LLNL under contract No. 
W-7405-Eng-48. GER is supported by the research grant PICT 
03-13291 (ANPCT). XB acknowledges financial support by the 
Spanish Ministerio de Education y Ciencia, under project ESP 
2003-00852. 


1 http://astroe2.gsfc.nasa.gov/docs/astroe/prop_tools/xrs_mat.html, with the most recent updates included (as of August 4, 2004) 
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Fig. 2.— Left: Ratio of data to fitted continuum for net 50 ks exposures, assuming inner disk inclination angles of 5°, 35° and 65°. Right: Detail of the ratio of 
data to fitted continuum for net 50 ks exposures, assuming inner disk inclination angles of 33°, 35° and 37° with the drop of the broad line going from left to right. 
Data have been grouped for fitting and presentation purposes. 
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Table 1 

Disk parameters, with 90% errors, for various effective exposure times 


Parameter 

input 

50 ks 

10 ks 

5 ks 

r 

1.80 

1.80±0.010 

1.80±0.015 

1.80 ±0.02 

£edge (keV) 

7.20 

7.19 ±0.05 

7.20±0.10 

7.12 ± 0.15 

^narrow (keV) 

6.415 

6.414±0.001 

6.414±0.001 

6.416 ±0.002 

^narrow (cV) 

30 

29.3 ±0.6 

29.3 ±1.2 

32.0 ±2.0 

Disk inclination 

35 

35.0 ±0.3 

35.2 ±0.6 

35.3 ±0.7 



